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ABSTRACT 

The RNA world hypothesis describes a stage in the 
early evolution of life in which RNA served as 
genome and as the only genome-encoded catalyst. 
To test whether RNA world organisms could have 
used cyclic trimetaphosphate as an energy source, 
we developed an in vitro selection strategy for 
isolating ribozymes that catalyze the 
triphosphorylation of RNA 5'-hydroxyl groups with 
trimetaphosphate. Several active sequences were 
isolated, and one ribozyme was analyzed in more 
detail. The ribozyme was truncated to 96 nt, while 
retaining full activity. It was converted to a trans- 
format and reacted with rates of 0.16min~^ under 
optimal conditions. The secondary structure 
appears to contain a four-helical junction motif. 
This study showed that ribozymes can use 
trimetaphosphate to triphosphorylate RNA 
5'-hydroxyl groups and suggested that RNA world 
organisms could have used trimetaphosphate as 
their energy source. 

INTRODUCTION 

The early evolution of life went through a stage in which 
RNA molecules served as genome and as genome-encoded 
catalyst, according to the RNA world hypothesis (1-3). 
This hypothesis is widely accepted due to several lines of 
evidence [for recent reviews, see (4,5)]. Direct fossil 
evidence of the RNA world, however, may never be 
found because the expected chemical components of these 
organisms would almost certainly not have survived more 
than three billion years until today (6). To determine how 
an RNA world organism could have functioned, re- 
searchers are attempting to re-create an RNA world 
organism in the laboratory (7-12). The replication of an 
RNA world organism would have relied on the synthesis 
of nucleotides (13-19), the generation of chemically 
activated nucleoside 5'-phosphates (see later in the text) 
and their polymerization (10). The focus of this study is 
the generation of chemically activated 5'-phosphates. 



The polymerization of RNA is entropically unfavorable 
in aqueous environments. Strategies to energetically drive 
RNA polymerization may have been different in two 
phases of the RNA world. First, in the prebiotic phase 
and perhaps in an early phase of the RNA world, RNA 
polymerization may have been driven by the evaporation 
of water from a solution (20) and from lipid encapsula- 
tions (21), or by the activation of the nucleotide's 5'-phos- 
phate with a wide range of leaving groups (22) such as 
adenine (23), cyanide (24), imidazole (25,26), or 2- 
methylimidazole (27). Importantly, these activation 
groups lead to hydrolysis in aqueous environment within 
hours. Second, in later phases of the RNA world, when 
lipid vesicles presumably encapsulated aqueous droplets of 
self-replicating RNA systems (9,28-30), kinetically more 
stable activation groups would have been important. This 
requirement is fulfilled by nucleoside triphosphates, the 
universal energy currency in all known hfe forms. 

Nucleoside triphosphates can be generated from nucleo- 
sides and cyclic trimetaphosphate (TMP) (31). TMP is a 
prebiotically plausible compound because it can be 
generated by volcanic activity (32), by the erosion of phos- 
phide minerals (33) and by heating of phosphate in the 
presence of urea (34). Additionally, TMP results from the 
self-reaction of the termini of hnear polyphosphates in 
aqueous medium (35). The prebiotic availability of TMP 
is supported by the recent finding of abundant, reactive 
and reduced phosphorus species in sediments of 3.5 
billion-year-old marine sediments (36). These phosphorus 
species can form polyphosphates including TMP if sources 
such as hghtening provide hydroxyl radicals (33). 
Although TMP is the most reactive polyphosphorylating 
reagent of alcohols among all polyphosphates (37), the 
reaction between nucleoside 5'-hydroxyl groups and 
TMP proceeds efficiently only at pH values >12 (31). At 
such high pH, RNA world organisms would hydrolyze 
quickly. Therefore, RNA world organisms would require 
a catalyst to use TMP as energy source. 

The triphosphorylation of nucleoside 5'-hydroxyl 
groups with TMP appears to be well within the range of 
catalysis by ribozymes. Most ribozymes are acid-base 
catalysts and have a particular capacity for phosphoryl 
transfer reactions (38-47). To identify ribozymes for new 
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reactions, the only known route is the in vitro selection 
from large pools of random RNA sequences (42,48,49). 
Therefore, we aimed to develop an in vitro selection 
scheme for the isolation of ribozymes that catalyze the 
triphosphorylation of RNA 5'-hydroxyl groups using 
TMP. 

Here, we describe the in vitro selection of a ribozyme 
that uses TMP to triphophorylate RNA 5'-hydroxyl 
groups. After five to eight rounds of selection, we 
obtained several active ribozymes, one of which was 
analyzed in more detail. This ribozyme was able to 
triphosphorylate RNAs in cis and in trans, with a kcat of 
up to 0.16min~'. These results suggested that RNA world 
organisms could have used TMP as their energy source. 

MATERIALS AND METHODS 

In vitro selection 

The DNA pool for the selection was generated from a 1 87- 
nt long DNA oligomer ('ultramer'. Integrated DNA 
Technologies (IDT)) that included 150 positions with 
randomized sequence (phosphoramidites were hand 
mixed) and primer binding sites on both termini. This 
187-mer was annealed to a 79-mer, which attached the 
sequence of the hammerhead ribozyme and the promoter 
of T7 RNA polymerase to the 5' -terminus of the pool. 
Amplification with short primers generated double- 
stranded copies of the pool. The DNA pool was transcribed 
into the RNA pool under standard conditions, during 
which the 5'-terminal hammerhead ribozymes cleaved 
themselves from the pool, generating 5'-hydroxyl groups 
on the RNA pool. The sequence of the transcript was 5'- 
GGGCGGTCTCCTGACGAGCTAAGCGAAACTGCG 
GAAACGCAGTC GAGACCGAGATGTT-N 1 50-CGCCA 
GTTAAGCTCCAGC-3', where the removed hammerhead 
ribozyme sequence is underlined. The RNA pool was 
purified by denaturing polyacrylamide gel electrophoresis 
to remove any uncleaved pool molecules. 

The 5'-hydroxylated RNA pool was incubated with TMP 
in a buffer containing 100 mM MgCli, 50 mM trisodium 
TMP (freshly dissolved and sterile filtered) and 50 mM 
Tris-HCl, pH 8.3. All triphosphorylation reactions during 
the selection were done at pH 8.3. When the solutions with 
TMP and MgCl2 were combined, the pH dropped to ~4.5. 
This pH was restored to the pH of the TMP solution alone 
(pH ~6) before the buffer (Tris-HCl, pH 8.3) was added. 
After 3h or 5min (depending on the selection round) of 
incubation at room temperature, RNAs were ethanol 
precipitated. The large salt pellet was extracted with a 
small volume of cold water to remove most of the salt. 
The remaining pellet was dissolved in water, desalted by 
size exclusion chromatography (P30 spin-columns; Bio- 
Rad), ethanol precipitated and redissolved in water. 

The recovered RNA pool molecules were heat 
renatured (2min/80°C) with a 1.25-fold molar excess of 
the R3C hgase ribozyme (50) (whose arms were designed 
complementary to the RNA pool 5'-terminus and the 
biotinylated capture oligonucleotide) and a 1.5-fold 
molar excess of biotinylated capture oligonucleotide 
[5'-biotin-d(GAACTGAAGTGTATG)rU-3'], in 100 mM 



KCl and 100 mM Tris-HCl, pH 8.5. The solution was 
diluted to 400 nM pool RNA, 500 nM hgase ribozyme, 
600 nM capture ohgonucleotide and adjusted to 50 mM 
KCl, 25 mM MgCb, 2mM spermidine, 20% (w/v) PEG 
8000 and 50 mM Tris-HCl, pH 8.5. After incubation for 
3h at 30° C, magnesium was chelated by an excess of 
Na2EDTA, and the mixture was heated (10min/50°C) 
with a 10-fold excess of a DNA complementary to the 
ligase ribozyme to free the ligated RNAs. The biotinylated 
nucleic acids were captured during 30min of agitation 
with streptavidin-coated magnetic beads (Promega) con- 
taining a 1.5-fold excess of biotin binding sites over 
biotinylated capture RNAs. Captured RNAs were 
washed first with 50 mM KCl, 20mM 4-(2- 
hydroxyethyl)-l-piperazineethanesulfonic acid (HEPES)/ 
KOH pH 7.2 and 0.01% Triton X-100, then with 
20 mM NaOH and 0.01% Triton X-100. RNA pool mol- 
ecules were eluted from the magnetic beads by heating 
(3min/65°C) with 95% formamide/1 mM NasEDTA. 
The effective complexity of the RNA pool was 1.7 x lO'"^ 
sequences, based on an initial complexity of the double- 
stranded DNA hbrary of 2.4 x lO' sequences, a total of 
1.6nmol of RNA pool molecules that entered the hgation 
step, losses of ~60% in the ligation step, losses of ~50% 
in the capture on streptavidin-coated beads and losses of 
~20% in further processing steps (data not shown). After 
ethanol precipitation, the RNAs were reverse transcribed 
using Superscript III reverse transcriptase (Invitrogen) 
and polymerase chain reaction (PCR) amplified using 
Taq polymerase. The 5'-PCR primer added a selective 
step because a part of its binding site was generated by 
the biotinylated capture ohgonucleotide. The DNA pool 
for the next round of selection was regenerated by PCR 
amplification with the 79-nt primer that was used to attach 
the hammerhead ribozyme and the T7 RNA polymerase 
promoter. We used a previously published protocol for 
mutagenic PCR (51). 

In pilot PCR amplifications of the reverse-transcribed 
products, short PCR artifacts created the dominant 
product after ~30 PCR cycles. We assumed that this 
was because some pool molecules allowed the PCR 
primers to anneal within the randomized region, and 
shorter amplicons are known to outcompete longer 
amplicons during PCR (52). To avoid that problem we 
'doped' the RNA pool molecules after the incubation 
with TMP with 1/37 000 molar amount of pool 0 mol- 
ecules that carried a 5'-triphosphate. Because these 
5'-triphosphorylated molecules were hgated efficiently to 
the capture oligonucleotides, the doped pools required 
only 18-30 PCR amplification cycles and suppressed 
short PCR artifacts. Although this may have slowed the 
selection compared with a protocol without PCR 
problems, the 1/37 000-fold did not seem to harm the en- 
richment over multiple rounds of selection. 

Triphosphorylation assays 

For the self-triphosphorylation assay based on the R3C 
hgase ribozyme, cloned ribozymes were first transcribed 
in vitro (53) and purified by denaturing polyacrylamide 
gel electrophoresis (PAGE). Triphosphorylation reactions 
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with 5 |iM triphosphorylation ribozyme were incubated 
with 50 mM TMP, lOOmM MgCl2 (corresponding to 
50 mM free Mg^^) and 50 mM Tris-HCl (final pH 8.1) 
for 3h at 22°C. The reaction mixture was diluted 10-fold 
with a buffer such that buffer concentrations were 500 nM 
ribozyme, 500 nM R3C ligase ribozyme, 500 nM of an 
ohgonucleotide (5'-[^^P]-d(GAACTGAAGTGTATG)rU- 
3'), 100 niM KCl, 100 mM Tris-HCl, pH 8.5, and 
15mM Na2EDTA. After heat denaturation (2min/ 
65°C), the mixtures were diluted another 2-fold and 
adjusted to 20% (w/v) PEG 8000, 22 mM free Mg^^, 
50 mM KCl, 2.5 mM spermidine and 100 mM Tris-HCl 
pH 8.5. This ligation reaction was incubated at 30° C for 
2h, stopped by ethanol precipitation and the products 
were separated by denaturing 10% PAGE. Gels were 
exposed to phosphorimager screens, scanned on a PMI 
phosphorimager (Bio-Rad) and signals were quantitated 
using the software quantity one. 

The self-triphosphorylation assay based on the 
DNAzyme used triphosphorylation ribozymes that were 
internally [''^P] labeled during transcription and purified 
by denaturing PAGE. Triphosphorylation reactions with 
8 |iM ribozyme were incubated as described earlier in the 
text. At specific time points, 2 |il aliquots of the reactions 
were quenched by adding them to 8 |.il of a solution con- 
taining 31.25 mM NaaEDTA, 375 mM NaCl, 32pmol of a 
DNAzyme that recognized the 14 nt conserved at all 
triphosphorylation ribozyme 5'-termini, and 32pmol of a 
DNA oligonucleotide that was complementary to the first 
25 nt of the 'Niso- sequence to improve accessibihty for the 
DNAzyme. After heat renaturation (2min/80°C), the 
DNAzyme reaction was initiated by adding 10|il of a 
solution containing 200 mM MgClz, lOOmM Tris-HCl 
pH 8.3, and incubated for 1 h at 37°C. Two microliter of 
these reactions were quenched with 8 |.il of PAGE loading 
buffer containing 16pmol of a DNA ohgonucleotide com- 
plementary to the DNAzyme, 90% (v/v) formamide and 
37.5 mM Na2EDTA. Products were separated by 
denaturing 22.5% PAGE and analyzed as described 
earlier in the text. 

The ?ra«.v-triphosphorylation assay used substrate mol- 
ecules that were internally [^^P] labeled during transcrip- 
tion. Two hammerhead ribozymes were contained in the 
transcript, one at the 5'-terminus and one at the 
3'-terminus. The 5'-terminal hammerhead ribozyme 
generated the 5'-hydroxyl group, whereas the 3'-terminal 
hammerhead ribozyme generated a homogeneous 
3'-terminus (54). Cleaved transcripts were purified by 
denaturing PAGE. Reaction mixtures with 5.5 |iM 
triphosphorylation ribozyme and substoichiometric con- 
centrations of radiolabeled substrate were incubated as 
described earher in the text. At specific time points, 
1.5 yil of the reactions were added to 6.5 |il of formamide 
PAGE loading buffer containing 20 mM Na2EDTA, heat 
renatured and separated by denaturing 20% PAGE. 
Quantitation of the signals was as aforementioned. For 
the titration of the TMP concentration, the chelation of 
Mg2+ with TMP generated a sHght reduction in the pH for 
the reaction with 200 mM [TMP] such that the pH of 8.1 
(all other reactions) was reduced to pH 7.6. 



TMP was obtained as trisodium salt from MP 
Biomedicals (cat # 150160; lot # 9425H). The supplier 
fisted 0.4 mg/kg of arsenic and <4 mg/kg of lead as contam- 
ination. We tested for other polyphosphates in the prepar- 
ation by electrospray-ionization Mass Spectroscopy 
(ESI-MS) analysis. The results suggested that the used 
TMP sample contained -80% TMP, 10% finear triphos- 
phate and 10% of finear tetraphosphate, under the assump- 
tion that these molecules have similar ionization efficiencies 
under ESI-MS. Traces of orthophosphate and no detect- 
able amounts of pyrophosphate or polyphosphates longer 
than tetraphosphate were found. 

Generation of modified ribozyme constructs 

Sequence modifications were introduced into ribozyme 
constructs by PCR mutagenesis or the Quickchange 
protocol for site-directed mutagenesis (Stratagene). AU 
DNA sequences were confirmed by cloning into pUC19 
and sequencing. The only exceptions were ribozyme con- 
structs where mutations near the 3'-terminus could be 
introduced during PCR amplification of the templates for 
transcription. RNAs were generated by runoff transcrip- 
tion by T7 RNA polymerase as described (53) and purified 
by denaturing PAGE. The sequences in the mutated 
substrate recognition duplex of the /ra«.y-reaction were 
5'-GAGACCGUACAUUU-3' for the substrate and 5'- 
GGAUGUAC ... -3' for the ribozyme 5'-terminus. 
Underfined nucleotides differ from the wild-type sequence. 

Mass spectrometric analysis 

Substrate RNAs were prepared as described for the trans- 
triphosphorylation reaction, but without internal 
radiolabefing. A fraction of the substrate was incubated 
with an excess of the TPRl ribozyme for 3 h. The product 
of the triphosphorylation reaction and unreacted substrate 
RNA were purified by denaturing PAGE, and desalted on 
Ci8 Zip-tips (MiUipore). About 20pmol of each sample 
were dissolved in 4|.il containing 3-hydroxypicohnic acid 
and dianimonium hydrogen citrate, and spotted on 
matrix-assisted laser desorption/ionization (MALDI) 
targets. The mass spectra were recorded by Dr 
Yongxuan Su at the UCSD Molecular Mass 
Spectrometry Facility in negative ion mode on a Bruker 
Biflex IV, MALDl-time-of-flight mass spectroscopy 
(TOFMS). 

SHAPE analysis 

For SHAPE analysis, the TPRl ribozyme was extended at 
its 3'-terminus by 12nt such that a 12-nt reverse transcrip- 
tion primer could facifitate analysis of the complete 
ribozyme. The fra«.v-reacting TPRl ribozyme was 
incubated at a concentration of 1 \\M with a 1.5-fold stoi- 
chiometric excess of its substrate, in lOOmM MgCl2, 
50 mM trisodium TMP and 50 mM HEPES/KOH pH 
8.0. To 9 |il of this solution, 1 |.il of dimethyl sulfoxide 
(DMSO) or a 20 mM solution of 1M7 (55) in DMSO 
was added. After incubation for 3min at room tempera- 
ture, the reaction products were ethanol precipitated. 
Reverse transcription with a 5'-[^^P]-radiolabeled 12-nt 
long reverse transcription primer was facilitated by 
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Superscript III reverse transcriptase (Invitrogen). 
Reaction products were ethanol precipitated, redissolved 
in formamide PAGE loading buffer, heated (2min/80°C) 
and separated by denaturing 10% and 20% PAGE. Signal 
quantitation was done as described for the 
triphosphorylation assays. Only primary reaction sites 
were observed because >80% of the reverse transcription 
products had full length, and signals appeared with consist- 
ent kinetics in a concentration series of 1M7 (not shown). 
The signals were judged as strong, or weak when their in- 
tensity was at least 30%, or 8.5% of the strongest SHAPE 
signal intensity, respectively. Markers were synthetic DNA 
molecules with 5'-[ '^P] label that had the exact sequence of 
expected reverse transcription products and a length of 25, 
47 and 70 nt. 

RESULTS 

In vitro selection 

An in vitro selection procedure was established to obtain 
ribozymes that use TMP to triphosphorylate RNA 
5'-hydroxyl groups (Figure 1). First, a double-stranded 
DNA pool was generated that encoded the T7 RNA poly- 
merase promoter, a hammerhead ribozyme and 150 
randomized nucleotides flanked by two constant primer 



binding sites. The hammerhead ribozyme served to 
generate 5'-hydroxyl groups at the 5'-termini of the pool 
molecules (56). The hammerhead ribozyme cleaved itself 
efficiently off the pool under transcription conditions so 
that >95% of the pool was cleaved after 2h of transcrip- 
tion incubation (data not shown). The gel-purified pool 
molecules were then incubated for 3h with 100 mM 
MgCl2, 50 mM TMP and 50 mM Tris-HCl pH 8.3 to 
allow catalytically active pool molecules to 
triphosphorylate their own 5'-terminus. We found that 
TMP chelated an equimolar amount of Mg'^^ (data not 
shown), therefore the triphosphorylation conditions 
provided 50 mM free magnesium ions. To isolate pool 
molecules that carried a 5'-triphosphate, the R3C ligase 
ribozyme was used (50). This ribozyme created 3'-5'- 
phosphodiester bonds between 5'-triphosphorylated pool 
molecules and the 3'-hydroxyl group of a 5'-biotinylated 
oHgonucleotide. The biotin modification allowed 
capturing the pool on streptavidin-coated magnetic 
beads, and stringent washing of pool molecules that 
were coupled covalently to the biotin moiety. The 
captured pool molecules were reverse transcribed and 
PCR amplified. In a second PGR step, the T7 RNA poly- 
merase promoter and the 5'-terminal hammerhead 
ribozyme sequence were regenerated, completing one 
round of the selection. In a mock selection round, we 
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Figure 1. Scheme for the in vitro selection of triphosphorylation ribozymes. (A) A double-stranded DNA library containing the promoter for T7 
RNA polymerase (T7), the sequence for a hammerhead ribozyme (HhRz) and a randomized sequence (N150) was transcribed into RNA. (B) The 
5'-terminal hammerhead ribozyme cleaved itself off cotranscriptionally (filled triangle), generating a 5'-terminal hydroxyl group on the RNA library. 

(C) The RNA library was incubated in the presence of TMP such that active ribozymes could triphosphorylate their 5'-terminus. 

(D) Triphosphorylated RNA molecules were reacted with the 3'-hydroxyl group of a short biotinylated RNA, by the R3C ligase ribozyme. 

(E) Ligated RNAs were captured via their biotin modification on streptavidin-coated magnetic beads and washed stringently. The RNAs were 
then (F) reverse transcribed and (G) amplified by PCR to generate the DNA pool for the next round of selection. 
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estimated that 5'-triphosphorylated RNAs were enriched 
at least 9000-fold in a single selection round 
(Supplementary Figure SI). 

Based on the effective complexity of the RNA pool of 
~1.7xl0''* sequences (see 'Materials and Methods' 
section) and the enrichment factor of ~10'^-fold/selection 
round (Supplementary Figure SI), we anticipated that 
catalytically active sequences could dominate the RNA 
pool after four selection rounds. To follow the progress 
of the selection, we monitored the number of PCR cycles 
that were required after each reverse transcription reaction 
to generate clearly detectable bands on agarose gels 
(Supplementary Figure S2). Satisfyingly, the number of 
PCR cycles dropped from at least 17 cycles in rounds 
one to 3-8 cycles in round four. Two modifications were 
used for the later selection rounds. First, mutagenic PCR 
was introduced in the amplification steps to allow the pool 
molecules to explore their sequence neighbors and perhaps 
find more active sequences. Second, the population was 
spht into two branches, one exposed to the same selection 
pressure as before (3 h of incubation with TMP), and one 
exposed to high selection pressure (only 5 min of incuba- 
tion with TMP). In round seven, a crossover branch was 
added to allow clones from low selection pressure to cross 
into the high selection pressure branch of round eight. 
A total of 36 sequences were obtained from the pools 
after five and eight selection rounds, from low and high 
selection pressure branches (Supplementary Figure S3). 
The clones were named after their selection history: R5 
or R8 denoting their round of isolation, a combination 
of the numbers 3 and 5 denoting the incubation times 
with TMP (3h and/or 5 min) and C with the clone 
number. 



Identification of the most active ribozyme clones 

To identify ribozynies with high activity, we screened the 
36 isolated clones for activity. This screen used the same 
procedure as during the selection, reacting the ribozyme 
with TMP followed by hgation to a short RNA using the 
R3C hgase ribozyme, but a 5'-radiolabeled ohgonucleo- 
tide was used instead of a biotinylated oligonucleotide 
(Supplementary Figure S4). The fraction of shifted oligo- 
nucleotides in a denaturing gel reported the hgation effi- 
ciency and indirectly the triphosphorylation efficiency. 
The 36 ribozyme clones were sorted according to their 
activity, which revealed a wide distribution of activities 
(Figure 2). Several sequences appeared multiple times, 
but their multiplicity was not correlated with increased 
activity. The activity of several ribozyme clones in this 
assay was ~2-fold higher than that of a positive control, 
made from pool RNA containing a 5'-triphosphate. We 
assumed that this was because the assay was not directly 
monitoring triphosphorylation activity, and the clones 
with peculiarly high activity worked especially well 
together with the R3C hgase ribozyme used in the selec- 
tion procedure. To identify those ribozyme clones that 
were most efficient in catalyzing the triphosphorylation 
reaction, we chose the eight most active clones from the 
hgase ribozyme-based assay. 

To directly monitor the triphosphorylation status at the 
ribozyme 5'-termini, a different gel shift assay was used 
(Figure 3). Because the gel shift of the triphosphate 
group was too small to detect migration differences of 
the 182-nt long ribozymes, we used the 8-17 DNAzyme 
(57) to cleave off the eight 5'-terminal nucleotides. The 
triphosphorylation ribozymes were internally labeled; 
therefore, the triphosphorylation status of the 8-nt 
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Figure 2. Screen of 36 in vitro selected RNAs for self-triphosphorylation activity, using an assay based on the R3C ligase ribozyme (see 
Supplementary Figure S4). The column heights show the percent of 5'-["^^P]-radiolabeled oligonucleotide that was ligated to the individual 
selected RNAs by the ligase ribozyme, after the selected RNAs were incubated with TMP for 3 h to allow for self-triphosphorylation. Therefore, 
the percent of ligated oligonucleotide was an indirect measure of self-triphosphorylation activity. As positive control, the unselected RNA pool was 
ligated to the radiolabeled oligonucleotide, facilitated by a 5'-triphosphate that was incorporated during transcription (white column and horizontal 
dashed line). The labels on the x-axis denote the clone names (round of isolation_branch of the evolution_clone number). The clones were sorted 
according to the average percent of ligated oligonucleotide. The sequences of all 36 clones are shown in Supplementary Figure S3. Classes of related 
sequences are indicated by symbols, with einpty squares (class 1), empty triangles (class 2), empty circles (class 3). filled circles (class 4), filled squares 
(class 5) and filled diamonds (class 6) each belonging to one class. Class 5 is represented by a single clone because the second clone had the identical 
sequence. No other clone showed related sequences among the 36 clones. The eight clones with the highest activity (arrows) were chosen for further 
analysis. Error bars are standard deviations from three independent experiments. 
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Figure 3. Kinetic analysis of the eight most promising ribozyme clones. (A) Schematic of the assay using the 8-17 DNAzyme. After reaction with 
TMP, internally [^^"P]-radiolabeled ribozymes were cleaved by the DNAzyme. This freed the eight 5'-terminal nucleotides, facilitating gel separation 
of triphosphorylated and unreacted RNAs. (B) Autoradiogram of products after the DNAzyme reaction and separation by denaturing 22% PAGE. 
An RNA that was not exposed to TMP (5'-OH) and an RNA that was transcribed with a 5'-terminal triphosphate {5'-PPP) were used as negative 
and positive controls, respectively. The incubation times with TMP are indicated on the top. The long fragment of the cleaved ribozymes (174nt) and 
possible remaining uncleaved ribozymes (182nt) migrated much slower than the 8-nt fragments. The 8-nt fragments were separated based on their 
phosphorylation status. The particular autoradiogram was from analysis of clone R8_35C18A. (C) Determination of triphosphorylation kinetics 
from signals as shown in (B). The percent of triphosphorylation of the 8-mer was plotted as function of the incubation time with TMP. Single- 
exponential fits are shown as black or gray lines. Error bars are standard deviations from three experiments. Syinbols are explained in (D). Single- 
exponential curve fits are shown in black lines for filled symbols and gray lines for empty symbols. (D) Symbols and clone names used in (C), 
together with the parameters obtained by curve fits, the maximal percentage of reacted ribozyme (Max.) and the observed pseudo-first order rate 
constant (kobs). 



fragment could be followed after separation on 
denaturing polyacrylamide gels and phosphorimaging. 
The 8-nt fragments of unreacted ribozymes coniigrated 
with RNAs containing a 5'-hydroxyl terminus, and 
the fragment after incubation with TMP comigrated 
with RNAs containing a 5'-triphosphate. This confirmed 
that the products were triphosphorylated at their 
5'-terminus. Quantitation of the signals after different 
triphosphorylation reaction times resulted in single- 
exponential time dependence, reveahng pseudo-first 
order rate constants between 0.01 and 0.03 min~' for 
each of the eight analyzed ribozymes. Although most ribo- 
zymes reacted to a final extent of 30-50%, one ribozyme 
(clone R5_3C21) reacted to >80%. This clone was chosen 
for further analysis. 

If the ribozymes were able to triphosphorylate single 
free nucleosides, they could directly generate nucleoside 
triphosphates. To test whether the eight most active ribo- 
zymes could triphosphorylate free nucleosides, they were 



truncated at their 5'-terminus by one nucleotide (to accom- 
modate the free nucleoside) and incubated with C'"*- 
radiolabeled guanosine. The products were separated by 
Polyethylenimine (PEI) cellulose thin layer 
chromatography (TLC) analysis and analyzed by autoradi- 
ography (Supplementary Figure S5). None of the tested 
ribozymes generated a detectable signal for guanosine 
triphosphate (GTP). The result was the same when the 
ribozymes were not truncated, or when they were incubated 
at pH 9.5 or at different temperatures (40°C and -20°C) 
(data not shown). 

Truncation of ribozyme R5_3C21 

We removed fragments from different portions of the 
182-nt long R5_3C21 ribozyme to test whether a shorter 
fully active ribozyme could be obtained (Figure 4). The 
activity of the resulting truncated ribozymes was measured 
with the gel shift assay that directly monitored the 
5'-triphosphorylation with the help of the 8-17 
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Figure 4. Truncation analysis of the ribozyme clone R5_3C21. (A) Schematic of the tested truncations. The numbering is relative to the complete 
sequence of the 182-nt long initial isolate (Supplementary Figure S3). Dotted lines indicate internal segments that were removed. Symbols to the right 
of each construct are consistent with the symbols in (B). (B) Triphosphorylation kinetics for the initial isolate and five truncated sequences. The 
kinetics was determined with the DNAzyme assay (Figure 3). Black lines represent single-exponential curve fits to the data. The symbols of 
truncations 1 and 2 are overlapping near the x-axis. 



DNAzyme. The choice of deleted ribozyme fragments was 
based on tlie predicted secondary structures of the full- 
length ribozyme (58). When portions of the 5'-domain of 
the ribozyme were deleted, triphosphorylation activity was 
abohshed (truncation 1 and 2). When portions of the 3'- 
domain were deleted, full triphosphorylation activity was 
retained (truncations 3 and 4). This allowed truncating the 
ribozyme to 96 nt, whereas further truncation to 92 nt 
caused a shght reduction in activity (truncation 5). 
Therefore, the further experiments used truncation 4 of 
ribozyme R5_3C21, with a length of 96 nt. This 
ribozyme was termed TPRl. 

Triphosphorylation of short RNAs in trans 

We were interested to see whether the TPRl ribozyme 
would be able to triphosphorylate a short substrate 
RNA in trans (Figure 5). For this purpose, a short 
hairpin structure near the ribozyme 5'-terminus was 
separated at its loop, and the duplex was elongated to 
different lengths (6-9 bp). This was done using substrates 
with a length of 12, 13, 14 and 15 nt that showed comple- 
mentarity to the corresponding ribozyme 5'-termini. After 
incubating the radiolabeled substrates with the corres- 
ponding TPRl ribozymes, the products were separated 
on denaturing polyacrylamide gels and the signals 
quantitated. The results showed that the TPRl ribozyme 
was triphosphorylating the substrates in trans. The 
shortest substrate, which formed a duplex of 6 bp with 
the ribozyme, resulted in slow reaction kinetics. In 
contrast, longer substrates that formed duplexes of 
7-9 bp showed the same kinetics and reacted to the same 
extent of reaction as the cw-acting TPRl ribozyme. The 
duplex length of 8 bp was chosen for further experiments, 
corresponding to a 14-nt substrate RNA. When 6 of the 
8 bp in this duplex were modified (see 'Materials and 
Methods' section for the sequences), the reaction kinetics 
did not change significantly (0.023 min~' in the mutant 
compared with 0.020 min~' in the wild-type). This 
suggests that the existence of the duplex but not the 
identity of its paired bases was important to mediate the 
trans-ve-dcXion. 



Several lines of evidence suggested that the TPRl 
ribozyme was generating 5'-triphosphate groups: the 
reaction between nucleoside 5'-hydroxyl groups and 
TMP was established previously (31), the reacted ribo- 
zymes could be hgated by the R3C ligase ribozyme 
(Supplementary Figure S4) and the ribozyme reaction 
products showed the same gel shift as a 5' -triphosphate 
(Figure 3B). In addition, the /ra«.v-reaction presented a 
good opportunity for a more rigorous test. The mass of 
the short substrate RNA was analyzed before and after 
the reaction with ribozyme and TMP (Figure 5D). The 
expected mass difference was 239.94 Da. Satisfyingly, 
the observed mass difference (239.83 Da) was in excellent 
agreement with the prediction, thereby confirming that the 
in vitro selected ribozymes used TMP to triphosphorylate 
the RNA 5'-hydroxyl groups. 

Secondary structure analysis 

To identify single-stranded nucleotides in the ribozyme, 
we subjected it to SHAPE analysis with l-methyl-7- 
nitroisatoic anhydride (1M7) (59). This reagent covalently 
reacts with 2'-hydroxyl groups that can adopt a specific 
geometry, and can therefore discriminate flexible nucleo- 
tides from double-stranded nucleotides, which is used to 
determine RNA secondary structures (60). The 1M7 
reactivities with our ribozyme identified several flexible 
regions and several regions that were probably base 
paired (Figure 6A). These constraints were used to build 
a secondary structure model that contained a central four- 
helix junction (Figure 6B). Four-helix junctions were 
identified previously for RNA molecules such as between 
the spliceosomal RNAs U2 and U6 (61) and in the 
internal ribosomal entry site of HCV (62). 

To test whether the two short helices formed, the pre- 
dicted base pairs between CIO and G33 on the ribozyme 
and between C5 on the substrate and G69 on the ribozyme 
were chosen for analysis (Figure 6C). At each base pair, 
three mutants were constructed and tested for activity. 
The single mutations removed the expected base pair 
and were expected to decrease activity. The double 
mutants were expected to restore the base pair and 
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catalytic activity. As predicted, there was no detectable 
activity for any of the single inutations ribozyine 
variants tested. The C5G/G69C double mutation 
restored -50% of the activity, but the C10G/G33C 
double mutation showed only sHght activity. This con- 
firmed the formation of the base pair between C5 and 
G69. It is possible that the C10/G33 base pair formed 
but that the identity of the bases was also important for 
activity. 

The sequence 5'-GGG-3' (G32-G34) could base pair to 
the sequence 5'-CCCCC-3' (C8-C12) in three possible 
registers, not only in the register described above (CIO/ 
G33). To test whether the hehx formed in any of the 
two alternative registers, base covariation experiments 
were performed at C9/G33 and C11/G33. No single 
mutant and no double mutant showed detectable activity 
(data not shown), in contrast to the register C10/G33. This 
suggested that the base pair between CIO and G33 was 
formed, and that the identity of CIO and/or G33 was also 
required for maximum activity. 

Dependence of kobs on |TMP|, |Mg2+| and pH 

To identify the reaction conditions at which the TPRl 
ribozyme worked most efficiently, the dependence of the 
observed reaction rate on the concentration of TMP, free 
magnesium ions and pH were investigated (Figure 7). 
TMP appeared to chelate an equimolar amount of Mg^^ 
(data not shown); therefore the excess of [Mg^"^] over 
[TMP] in our experiments was called the free [Mg^^]. 



When the concentration of free magnesium ions was 
kept constant at 50 mM, an increase in the TMP concen- 
tration led to an increase in the observed rate to 
0.03iTiin"' at lOOinM [TMP], at pH 8.1. The [TMP] de- 
pendence was consistent with a 1 : 1 stoichiometry between 
ribozyme and TMP (curve fit in Figure 7A). Increases in 
the magnesiuin concentration at a constant TMP concen- 
tration of 100 mM led to a further increase in the reaction 
rate, with a kobs of 0. 1 6 min" ' at 400 mM of free [Mg^^j, at 
pH 8.1 (Figure 7B). The pH dependence of the reaction 
was tested under these conditions [lOOmM [TMP] and 
400 mM free [Mg^^]], as well as under the reaction con- 
ditions used during the selection [50 mM [TMP] and 
50 mM free [Mg^"^]] (Figure 7C). In these reactions 
three buffers were used: 2-(N-morpholino)ethanesulfonic 
acid (MES)/NaOH, HEPES/NaOH and Tris-HCl. For 
all reactions, the observed triphosphorylation rate 
increased Hnearly with the pH in the region between 
pH 5.5 and 8.5, with a slope of about one. This is con- 
sistent with a single deprotonation step being limiting for 
the reaction, perhaps the deprotonation of the RNA 
5'-hydroxyl group in preparation for the nucleophilic 
attack on TMP. 



DISCUSSION 

We established an in vitro selection strategy for ribozymes 
with a novel activity, the triphosphorylation of RNA 
5'-hydroxyl groups. The procedure identified more than 
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Figure 6. Secondary structure analysis for the /fnH.s-triphosphorylating TPRl ribozyme. (A) Autoradiograms of SHAPE analysis products, after 
reacting the ribozymes with 1M7 and subjecting them to reverse transcription with a 5'-[^^"P]-radiolabeled primer. The left image shows a separation 
by 20% PAGE to separate short products, whereas the right image shows a separation by 10% PAGE to show the longer products. Each image 
shows three lanes, where (— ) denotes a negative control with DMSO, (+) denotes the SHAPE reaction with 1M7 dissolved in DMSO and M denotes 
a marker lane with three 5'-radiolabeled DNAs that have the identical sequence as the expected reverse transcription products. The primer has a 
length of I2nt, base paired to 12nt that were appended to the ribozyme 3'-terminus. The SHAPE signal is shifted by one nucleotide relative to the 
length of the reverse transcription product because the reverse transcriptase stops at the nucleotide before the SHAPE modification. (B) The 
secondary structure was based on two types of analysis. Filled triangles indicate a strong signal from SHAPE analysis, whereas empty triangles 
indicate weaker signals. (C) Two base pairs suggested by the SHAPE data were tested by single mutations of each base partner, and double mutation 
that should have restored activity for correct base pairs. The observed reaction rate was determined as in Figure 5 and given for each mutant. 



a dozen independent ribozymes. One ribozyme was 
characterized in more detail. This ribozyme was able to 
work in trans, and its product was confirmed by mass 
spectroscopy. Its secondary structure appears to include 
a four-way helical junction. Its optimal reaction rate was 
at lOOmM TMP, 400 mM free [Mg^^] and increased pH. 

The observed rate with the TPRl ribozyme was 
0.16min~' under optimal conditions (100 mM [TMP], 
400 mM free [Mg'^^], pH 8.1) (Figure 7). As comparison, 
the rate of the uncatalyzed reaction (kuncat) in the absence 
of Mg^"^ was estimated to be 3.4 x 10~'*min~', at 500 mM 
[TMP] and pH 12 [Figure 2 in (31)]. Because the pH is 
limiting for the uncatalyzed reaction, its rate at pH 8.1 was 
assumed to be lO^'^-fold slower, at 4.3 x 10~ min~'. In 
addition, higher concentrations of TMP lead to faster re- 
actions for the uncatalyzed reaction (31). Because the in- 
fluence of [TMP] on kuncat may be non-hnear, we assumed 
a 3-fold lower k^mcat at lOOmM [TMP] compared with 



500mM [TMP]. This rate of 1.4 x lO'^min"' is 
1.1 X lO'^-fold below the rate of 0.16 min for the 
ribozyme-catalyzed reaction. Mg^"^ ions accelerate the 
uncatalyzed reaction (63) but no kinetic data are available 
for this Mg^"^ acceleration. Therefore, the value for k^at/ 
kuncat of ~10^ describes the rate enhancement of the 
ribozyme-catalyzed reaction in the presence of magnesium 
relative to the rate of the uncatalyzed reaction without 
magnesium. 

None of the six tested ribozymes showed 
triphosphorylation activity of free nucleosides 
(Supplementary Figure S5). It may be possible to 
overcome this in future experiments because many more 
ribozyme clones from our selection can now be tested for 
this activity. The fact that most of the isolated 36 clones 
appeared as single clones (Supplementary Figure S3) 
suggests that many more active sequences can be 
isolated. However, RNA world organisms could have 
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Figure 7. Dependence of reaction kinetics on reaction conditions, for the frajis-splicing TPRl ribozyme. (A) Dependence of the observed reaction 
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relied on the triphosphorylation of RNAs instead of free 
nucleosides to drive RNA polymerization: activation 
groups at the nucleoside 5'-phosphate are only necessary 
for RNA polyinerization in the 5'-3' direction. In 
contrast, RNA polymerization in the 3'-5' direction does 
not require nucleotide activation because chain elongation 
proceeds via the reaction of the activated RNA primer 
5'-terminus with the nucleoside 3'-hydroxyl group. This 
principle has been shown with an evolved variant of the 
class I ligase ribozyme that uses nucleoside triphosphates 
(NTPs) for the extension of RNA primers in both the 3'-5' 
and 5'-3' direction (64). This suggests that instead of NTP, 
a nucleoside could be used to extend the growing priiner 
5'-terminus. In the second step, the new RNA 5'-terminus 
could be triphosphorylated. Therefore, RNA polymeriza- 
tion in an RNA world organism could have proceeded via 
the alternating triphosphorylation of the RNA primer, 
and the extension of the primer 5'-terminus with a 
templated nucleoside. This polymerization in the 3'-5' dir- 
ection may have had an advantage for RNA world organ- 
isms: Nucleoside triphosphates are necessary for 5'-3' 
elongation, but their negative charge destabilizes the 
binding of the monoiner to the teinplate strand due to 
charge repulsion with the phosphodiester backbone of 
the template strand. This difference in affinity can be 
inferred froin measurements of stacking interactions in 
nucleosides and NTPs (65). In contrast, RNA polymeriza- 
tion in the 3'-5' direction could use the uncharged nucleo- 
sides for the primer extension step, with a corresponding 
increase in inonomer affinity to the template strand. 
Therefore, the alternating RNA 5'-triphosphorylation 
and extension with a nucleoside could have been an effi- 
cient RNA polymerization strategy for RNA world 
organisms. 

In addition to promoting RNA polymerization, RNA 
world organisms could have used RNA 5'-triphosphates 
for several other reactions. This could have included 
RNA-RNA ligation (42,50,66,67), the formation of 



phosphoamidate bonds (45), RNA capping (46,68) and 
ainino acid activation (69). Together with the results of 
this study, this suggests that RNA world organisins could 
have used TJMP to thermodynamically drive a wide range 
of reactions. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 
ACKNOWLEDGEMENTS 

The authors thank Chengguo Yao for pilot experiments 
with the hammerhead ribozyme, Manny Ares for the 
generous gift of 1M7 for SHAPE analysis and Ram 
Krishnamurthy for helpful discussions. 

FUNDING 

National Institutes of Health [molecular biophysics 
training grant T32 GM008326 to J.E.IVI.]; Hellman 
Family foundation [awards 2011/12 and 2012/13 to 
U.F.M.]; This material is based upon work supported by 
the National Aeronautics and Space Administration 
under Grant/cooperative agreement [No. NNX13AJ09G] 
issued through the Science Mission Directorate, Research 
Opportunities in Space and Earth Sciences -201 1 (ROSES- 
2011), in Astrobiology/Exobiology (to U.F.M.). Funding 
for open access charge: NASA grant [NNX13AJ09G]. 

Conflict of interest statement. None declared. 
REFERENCES 

1. Woese,C.R. (1967) The Genetic Code. The Molecular Basu for 
Genetic Expres.sion. Harper & Row, New York. 

2. Crick,F.H.C. (1968) The origin of the genetic code. /. Mol. Biol, 
38, 367-379. 



Nucleic Acids Research, 2014, Vol. 42, No. 7 4777 



3. Orgel,L.E. (1968) Evolution of the genetic apparatus. /. Mol. 
Biol. 38, 381-393. 

4. Robertson.M.P. and Joyce.G.F. (2012) The origins of the RNA 
world. Cold Spring Harh. Perspect. Biol, 4, pii: a003608. 

5. SzostakJ.W. (2012) The eightfold path to non-enzymatic RNA 
replication. /. Syst. Chem., 3, 2. 

6. Wacey,D., Kilburn,M.R., Saunders,M., CHffJ. and 
Brasier,M.D. (2011) Microfossils of sulphur-metabolizing cells 
in 3.4-billion-year-old rocks of Western Australia. Nat. Geosci., 
4, 698-702. 

7. Von Kiedrowski.G. (1986) A self-replicating hexadeoxynucleotide. 
AngL'w. Chem. Int. Ed. Engl., 25, 932-935. 

8. Orgel,L.E. (2004) Prebiotic chemistry and the origin of the RNA 
world. Crit. Rev. Biochem. Mol. Biol, 39, 99-123. 

9. Hargreaves,W.R. and Deamer.D.W. (1978) Liposomes from ionic, 
single-chain amphiphiles. Biochemistry, 17, 3759-3768. 

10. Johnston,W.K., Unrau,P.J., Lawrence,M.S., Glasner,M.E. and 
Bartel,D.P. (2001) RNA-catalyzed RNA polymerization: accurate 
and general RNA-templated primer extension. Seience, 292, 
1319-1325. 

11. Szostak.J.W. (2011) An optimal degree of physical and chemical 
heterogeneity for the origin of life? Philo.i. Tran.s. R. Soc. Lond. B 
Biol. Sci., 366, 2894-2901. 

12. Sczepanski,J.T. and Joyce.G.F. (2012) Synthetic evolving systems 
that implement a user-specified genetic code of arbitrary design. 
Chem. Biol., 19, 1324-1332. 

13. Sanchez,R.A. and OrgekL.E. (1970) Studies in prebiotic synthesis. 
V. Synthesis and photoanomerization of pyrimidine nucleosides. 
/. Mot. Biol., 47, 531-543. 

14. Powner,M.W., Gerland,B. and Sutherland.J.D. (2009) Synthesis of 
activated pyrimidine ribonucleotides in prebiotically plausible 
conditions. Nature, 459, 239-242. 

15. Powner.M.W., Sutherland.J.D. and Szostak,J.W. (2010) 
Chemoselective multicomponent one-pot assembly of purine 
precursors in water. /. Am. Chem. Soc., 132, 16677-16688. 

16. Joyce,G.F., Schwartz.A.W., Miller,S.L. and Orgel,L.E. (1987) 
The case for an ancestral genetic system involving siinple 
analogues of the nucleotides. Proc. Nail Acad. Sci. USA, 84, 
4398^402. 

17. Krishnamurthy,R. (2012) Role of pK(a) of nucleobases in the 
origins of chemical evolution. Acc. Chem. Rc.i., 45, 2035-2044. 

18. Hud,N.V., Cafferty,B.J., Krishnamurthy.R. and WiUiams.L.D. 
(2013) The origin of RNA and "my grandfather's axe". Chem. 
Biol, 20, 466-474. 

19. Chen,M.C., Cafferty,B.J., Mamajanov,!., Gallego,!., Khanam,J., 
Krishnamurthy,R. and Hud,N.V. (2013) Spontaneous prebiotic 
formation of a beta-ribofuranoside that self-assembles with a 
complementary heterocycle. /. Am. Chem. Soc. (epub ahead of 
print, doi:10.1021/ja410124v). 

20. Verlander,M.S., Lohrmann.R. and Orgel,L.E. (1973) Catalysts 
for the self-polymerization of adenosine cyclic 2', 3'-phosphate. 
/. Mol. EvoL, 2, 303-316. 

21. Rajamani,S., Vlassov,A., Benner,S., Coombs,A., 01asagasti,F. and 
Deamer,D. (2008) Lipid-assisted synthesis of RNA-Uke polymers 
from mononucleotides. Orig. Life Evol. Biosph, 38, 57-74. 

22. Prabahar,K.J., Cole,T.D. and Ferris,J.P. (1994) Effect of 
phosphate activating group on oligonucleotide formation on 
montmoriUonite: the regioselective formation of 3',5'-linked 
ohgoadenylates. /. Am. Chem. Soc, 116, 10914-10920. 

23. Prabahar,K.J. and Ferris,J.P. (1997) Adenine derivatives as 
phosphate-activating groups for the regioselective formation of 
3',5'-linked ohgoadenylates on montmoriUonite: possible 
phosphate-activating groups for the prebiotic synthesis of RNA. 
/. Am. Chem. Soc, 119, 4330-4337. 

24. Lohrmann,R. and Orgel.L.E. (1968) Prebiotic synthesis: 
phosphorylation in aqueous solution. Science, 161, 64-66. 

25. Lohrmann.R. (1977) Formation of nucleoside 5'- 
phosphoramidates under potentially prebiological conditions. 
J. Mol. Evol., 10, 137-154. 

26. Sawai,H. and Orgel,L.E. (1975) Letter: Ohgonucleotide synthesis 
catalyzed by the Zn-2+ ion. /. Am. Chem. Soc, 97, 3532-3533. 

27. Inoue.T. and Orgel,L.E. (1982) Oligomerization of (guanosine 5'- 
pliosphor)-2-methylimidazolide on poly(C). An RNA polymerase 
model. J. Mol. Biol., 162, 201-217. 



28. HanczycM.M., Fujikawa,S.M. and Szostak.J.W. (2003) 
Experimental models of primitive cellular compartments: 
encapsulation, growth, and division. Science, 302, 618-622. 

29. Chen.LA., Salehi-Aslitiani,K. and Szostak,J.W. (2005) RNA 
catalysis in model protocell vesicles. /. Am. Chem. Soc, 127, 
13213-13219. 

30. Budin.L, Debnath,A. and Szostak.J.W. (2012) Concentration- 
driven growth of model protocell membranes. /. Am. Chem. Soc, 
134, 20812-20819. 

31. Etaix,E. and Orgel.L.E. (1978) Phosphorylation of nucleosides in 
aqueous solution using trimetaphosphate: Formation of 
nucleoside triphosphates. /. Carhohydr. Nucleosides Nucleotides, 5, 
91-110. 

32. Yamagata,Y., Watanabe,H., Saitoh,M. and Nainba.T. (1991) 
Volcanic production of polyphosphates and its relevance to 
prebiotic evolution. Nature, 352, 516-519. 

33. Pasek,M.A., Kee,T.P., Bryant,D.E., Pavlov,A.A. and Lunine.J.I. 
(2008) Production of potentially prebiotic condensed phosphates 
by phosphorus redox cheinistry. Angew Chem. Int. Ed. Engl., 47, 
7918-7920. 

34. Osterberg.R. and Orgel.L.E. (1972) Polyphosphate and 
trimetaphosphate formation under potentially prebiotic 
conditions. /. Mol. Evol., 1, 241-248. 

35. Thilo,E., Schulz,G. and Wichmann,E.-M. (1953) Die Konstitution 
des Grahamschen und Kurrolschen Salzes. Zeilschrift fiir 
Anorganische und Allgemeine Chemie, 272, 182-200. 

36. Pasek,M.A., Harnmeijer,J.P., Buick,R., Gull.M. and Atlas,Z. 
(2013) Evidence for reactive reduced phosphorus species in the 
early Archean ocean. Proc. Natl Acad. Sci. USA, 110, 
10089-10094. 

37. Feldmann,W. (1967) Das Trimetaphosphat als 
Triphosphorylierungsmittel fiir Alkohole und Kohlenhydrate in 
waBriger Losung. Seine Sonderstellung unter den kondensierten 
Phosphaten. Chem. Ber., 100, 3850-3860. 

38. Kruger,K., Grabowski,P.J., Zaug.A.J., Sands,J., Gottschling.D.E. 
and Cech,T.R. (1982) Self-splicing RNA: autoexcision and 
autocyclization of the ribosomal RNA intervening sequence of 
Tetrahymena. Cell, 31, 147-157. 

39. Guerrier-Takada.C, Gardiner,K., Marsh,T., Pace,N. and 
Altman,S. (1983) The RNA moiety of ribonuclease P is the 
catalytic subunit of the enzyme. Cell, 35, 849-857. 

40. Hutchins,C.J., Rathjen,P.D., Forster,A.C. and Symons,R.H. 
(1986) Self-cleavage of plus and minus RNA transcripts of 
avocado sunblotch viroid. Nucleic Acids Res., 14, 3627-3640. 

41. Buzayan,J.M., Gerlach.W.L. and Bruening.G. (1986) Satellite 
tobacco ringspot virus RNA: A subset of the RNA sequence is 
sufficient for autolytic processing. Proc. Natl Acad. Sci. USA, 83, 
8859-8862. 

42. Bartel.D.P. and Szostak,J.W. (1993) Isolation of new 
ribozymes from a large pool of random sequences. Science, 261, 
1411-1418. 

43. Lorsch.J.R. and Szostak.J.W. (1994) In vitro evolution of new 
ribozymes with polynucleotide kinase activity. Nature, 371, 31-36. 

44. Wiegand,T.W., Janssen,R.C. and Eaton,B.E. (1997) Selection of 
RNA amide synthases. Chem. Biol., 4, 675-683. 

45. BaskervillcS. and Bartel,D.P. (2002) A ribozyme that ligates 
RNA to protein. Proc. Natl Acad. Sci. USA, 99. 9154-9159. 

46. Huang.F. and Yarus.M. (1997) 5'-RNA self-capping from 
guanosine diphosphate. Biochemistry, 36. 6557-6563. 

47. Unrau.P.J. and Bartel.D.P. (1998) RNA-catalysed nucleotide 
synthesis. Nature, 395. 260-263. 

48. Tuerk.C. and Gold.L. (1990) Systematic evolution of ligands by 
exponential enrichment: RNA ligands to bacteriophage T4 DNA 
polymerase. Science, 249. 505-510. 

49. Elhngton.A.D. and Szostak.J.W. (1990) /;; vitro selection of RNA 
molecules that bind specific ligands. Nature, 346, 818-822. 

50. Rogers,J. and Joyce,G.F. (2001) The effect of cytidine on the 
structure and function of an RNA ligase ribozyine. RNA, 1, 
395-404. 

51. Cadwell,R.C. and Joyce.G.F. (1992) Randomization of genes by 
PCR mutagenesis. PCR Methods Appl., 2, 28-33. 

52. Bustin,S.A. (2000) Absolute quantification of mRNA using 
real-time reverse transcription polymerase chain reaction assays. 
/. Mol. Endocrinol., 25, 169-193. 



4778 Nucleic Acids Research, 2014, Vol. 42, No. 7 



53. MilliganJ.F., Groebe.D.R., Witherell.G.W. and Uhlenbeck.O.C. 
(1987) Oligoribonucleotide synthesis using T7 RNA polymerase 
and synthetic DNA templates. Nucleic Acids Res., 15, 8783-8798. 

54. Price,S.R., Ito,N., Oubridge,C., AvisJ.M. and Nagai.K. (1995) 
Crystallization of RNA-protein complexes. I. Methods for the 
large-scale preparation of RNA suitable for crystallographic 
studies. /. MoJ. Biol., 249, 398^08. 

55. Turner,R., Shefer,K. and Ares,M. Jr (2013) Safer one-pot 
synthesis of the 'SHAPE' reagent l-methyl-7-nitroisatoic 
anhydride (lm7). RNA, 19, 1857-1863. 

56. Uhlenbeck,O.C. (1987) A small catalytic oligoribonucleotide. 
Nature, 328, 596-600. 

57. Santoro.S.W. and Joyce,G.F. (1998) Mechanism and utility 
of an RNA-cleaving DNA enzyme. Biochemistrv, 37, 
13330-13342. 

58. Zuker,M. (2003) Mfold web server for nucleic acid folding and 
hybridization prediction. Nucleic Acids Res., 31, 3406-3415. 

59. Mortimer,S.A. and Weeks,K.M. (2007) A fast-acting reagent for 
accurate analysis of RNA secondary and tertiary structure by 
SHAPE chemistry. J. Am. Chem. Soc, 129, 4144-4145. 

60. Hajdin,C.E., Bellaousov.S., Huggins,W., Leonard,C.W., 
Mathews,D.H. and Weeks,K.M. (2013) Accurate SHAPE-directed 
RNA secondary structure modeling, including pseudoknots. Proc. 
Natl Acad. Sci. USA, 110, 5498-5503. 

61. Sashital,D.G., Cornilescu,G., McManus,C.J., Brow,D.A. and 
Butclier,S.E. (2004) U2-U6 RNA folding reveals a group II 



intron-like domain and a four-helix junction. Nat. Struct. Mol. 
Biol, 11, 1237-1242. 

62. Berry,K.E., Wagliray,S., Mortimer.S.A., Bai,Y. and Doudna,J.A. 
(2011) Crystal structure of the HCV IRES central domain reveals 
strategy for start-codon positioning. Structure, 19, 1456-1466. 

63. Clieng,C., Fan,C., Wan,R., Tong,C., Miao,Z., Chen,J. and 
Zhao,Y. (2002) Phosphorylation of adenosine with 
trimetaphosphate under simulated prebiotic conditions. Orig. Life 
Evol. Biosph, 32, 219-224. 

64. McGinness,K.E., Wright,M.C. and Joyce,G.F. (2002) Continuous 
ill vitro evolution of a ribozyme that catalyzes three successive 
nucleotidyl addition reactions. Chem. Biol., 9, 585-596. 

65. Sigel,H. and Griesser,R. (2005) Nucleoside 5'-triphosphates: self- 
association, acid-base, and metal ion-binding properties in 
solution. Chem. Soc. Rev., 34, 875-900. 

66. Ikawa,Y., Tsuda,K., Matsumura,S. and Iiioue,T. (2004) De novo 
synthesis and development of an RNA enzyme. Proc. Nat! Acad. 
Sci. USA. 101, 13750-13755. 

67. Robertson,M.P. and Ellington.A.D. (1999) /« vitro selection of an 
allosteric ribozyme that transduces analytes to amplicons. Nat. 
BiotechnoL, 17, 62-66. 

68. Zaher,H.S., Watkins,R.A. and Unrau,P.J. (2006) Two 
independently selected capping ribozymes share similar substrate 
requirements. RNA, 12, 1949-1958. 

69. Kumar,R.K. and Yarus,M. (2001) RNA-catalyzed amino acid 
activation. Biochemistry, 40, 6998-7004. 



